arXiv:1502.06596v2 [hep-th] 18 Apr 2015 


NORDITA-2015-24 


Fractional Skyrmions and their molecules 

Sven Bjarke Gudnason^’^’j^ and Muneto Nitta^’|^ 

^Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 
^Nordita, KTH Royal Institute of Technology and Stockholm University, 
Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden 
^Department of Physics, and Research and Education Center for Natural Sciences, 
Keio University, Hiyoshi f-El, Yokohama, Kanagawa 223-8521, Japan 

(Dated: April 21, 2015) 

Abstract 

We study a Skyrme-type model with a quadratic potential for a field with S'^ vacua. We consider 
two flavors of the model, the first is the Skyrme model and the second has a sixth-order derivative 
term instead of the Skyrme term; both with the added quadratic potential. The model contains 
molecules of half Skyrmions, each of them is a global (anti-)monopole with baryon number 1/2. 
We numerically construct solutions with baryon numbers one through six, and find stable solutions 
which look like beads on rings. We also construct a molecule with fractional Skyrmions having the 
baryon numbers 1/3-I-2/3, by adding a linear potential term. 
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I. INTRODUCTION 


Half a century ago Skyrme made a proposal [I] that Skyrmions, i.e. topological objects 
characterized by the third homotopy group tts, could describe nucleons in the pion effective 
field theory [2] if augmented by a higher-order derivative term. This term was then later 
denoted as the Skyrme term. The Skyrme term is needed for stabilizing the Skyrmions 
against shrinkage. Later nucleons are known to be bound states of quarks, described by 
QCD. However, QCD at low energies is strongly coupled and hard to tackle and hence the 
Skyrme model, with its small number of parameters, has remained an attractive model at 
low energies. In the limit of a large number of colors, the Skyrmion is exactly the baryon 
[3]. It has also recently been used in holographic QCD [U [5]. 

It has recently been found that half Skyrmions stably appear as constituents of a lattice at 
finite baryon density [B]. It is therefore important to understand the physical consequences 
of a half or non-integer topological charge of half (or fractional) Skyrmions, such as the 
statistics of them. It is, however, not easy to pick up an isolated half Skyrmion because it 
is not stable in isolation. 

The fractionality of topological charge, however, has been better understood in lower 
dimensions. As a lower-dimensional analog of Skyrmions, baby Skyrmions were proposed in 
an 0(3) sigma model with a fourth-order derivative term in 2-|-l dimensions [H [Hj, which 
topologically are lumps characterized by t ^2 instead of tts, i.e. the case for usual Skyrmions. 
2-1-1 dimensional Skyrmions often appear in various condensed-matter systems such as fer- 
romagnets and quantum-Hall systems. A baby-Skyrme model with an XY-type potential 
V = (which is also called easy-plane in ferromagnets) admits half baby Skyrmions 

[MI]. Each of them can be regarded as a global vortex covering the northern or southern 
hemisphere of the target space and consequently having a half 7r2 topological charge. Each 
of them has logarithmically divergent energy when (if) isolated, because it is a global vortex. 
If one separates the two by an infinite distance, the energy diverges logarithmically and thus 
they are confined to the form of a molecule. If a U(l) subgroup of the 0(3) symmetry in 
the 0(3) model is gauged, each constituent baby Skyrmion becomes a local vortex (here 
the fourth-order derivative term is not needed for stability), (still) carrying half a unit of 712 
charge [EHE]. In this case, if we choose the gauge coupling and the scalar coupling to be 
the same, the baby Skyrmions are BPS and can be embedded in a supersymmetric theory 
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M- An entire molecule can be separated from other molecules without any cost of energy. 
In this phase the two half Skyrmions can be separated by a finite distance at a finite cost of 
energy. In this sense, the half Skyrmion can be isolated. 

In this paper, we construct 3+1 dimensional half Skyrmions and their molecules in two 
Skyrme models with a potential term in the form of 1^ = where we use the notation of 

the 0(4) sigma model ua^x) (A = 1, 2, 3,4) with the constraint = 1- This potential 

is a potential for the would-be vacuum moduli and we will denote it a Heisenberg-type 
potential. The two models we consider in this paper are the conventional Skyrme model and a 
Skyrme-like model with the fourth-order derivative term replaced by a sixth-order derivative 
term. The latter is inspired by the BPS Skyrme model [16]. The sixth-order derivative term 
is the baryon-current density squared [TfiHT^ and this term alone with an adequate potential 
provides the basis of an integrable subsector of the model. For a short-term notation, we 
call them the 2+4 and 2+6 model, respectively. Each constituent of the molecule is a half 
Skyrmion carrying half a baryon number, i.e. the topological charge of tts. This turns out to 
be the case for both the 2+4 and the 2+6 models. The constituents are more separated for 
the 2+6 model than for the 2+4 model. This is in contrast with the conventional Skyrmions 
(i.e. without our potential) for which the configuration for B = 1 is spherically symmetric, 
that for B = 2 is toroidal, and those for B > 2 have energy distributions with some point 
symmetry [ 201121 ]. We take also the limit of the mass-parameters of our potential going to 
zero and recover the spherically symmetric Skyrmion for 5 = 1 and the transition is smooth 
as expected. We construct also higher baryon numbered solutions numerically; specifically 
B = 2,3, 4, 5, 6 and find that the lowest-energy states are found using an axially symmetric 
initial guess (basically a torus) and the relaxation method finds the solutions which look like 
beads on rings. We confirm that the configurations of the form of beads on rings are really 
the lowest-energy states, by checking different initial conditions; i.e. both those with an 
axial symmetric Ansatz and others with a rational map Ansatz having an appropriate point 
symmetry [2(11 [2T] . The rings of higher B are all stable against decaying into a sum of smaller 
B rings (for the B we have found explicitly here). This is a three-dimensional analog of baby 
Skyrmions as beads on a ring in a baby-Skyrme model with the XY-potential [IDIIII]. We 
also construct fractional Skyrmions with an arbitrary baryon number by using the potential 
V = m?{n^ “ c)^, for which each unit Skyrmion charge is split into two fractional Skyrmions 
with baryon number (1 + c)/2 and (1 — c)/2, respectively. We call them unequal fractional 
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molecules. We also note that fractional Skyrmions can be identified as global monopoles 
having divergent energy in an infinite system. 

If one gauges an SO(3) subgroup acting on ni,n 2 ,n 3 (a diagonal subgroup of the chiral 
symmetry), a monopole becomes the local’t Hooft-Polyakov monopole having finite energy 
[22H29j, similarly to a vortex as a half lump in 2+1 dimensions [la na [ig. Our case is a 
global analog of this case. 

This paper is organized as follows. In Sec. |TT| we present our two models and their 
numerical solutions in Sec. imiLe. of fractional-Skyrmion molecules. In Sec.|IV[ we construct 
an isolated fractional Skyrmion as a global monopole. In Sec. |V| numerical solutions with 


higher baryon numbers B = 2,3,4, 5,6 are presented. Finally, in Sec. |VI[ we construct a 
molecule where the weight of the two constituents of the molecule is altered, i.e. an unequal 
fractional molecule. Sec. m is devoted to a summary and discussions. Two investigations 
are delegated to the appendices. In App. |^we try different initial guesses for the numerical 
relaxation and find only metastable solutions. Finally, App. |B| shows how the molecule turns 
into a spherical Skyrmion by turning off the potential under study. 


II. A SKYRME-LIKE MODEL WITH A HEISENBERG-TYPE POTENTIAL 

We consider the SU( 2 ) principal chiral model with the addition of the Skyrme term and 
a sixth-order derivative term in d = 3 + 1 dimensions. In terms of the SU(2)-valued field 
U{x) G SU( 2 ), the Lagrangian which we are considering is given by 

£ = -|tr {d^U^d^U) + C 4£4 + cg/le - V{U), (1) 

where we use the mostly-positive metric and the higher-derivative terms are given by 

jC, = -^tv{[U^d,U,U^dMY), ( 2 ) 

A = [U^dMU^d,UU^dM]Y , (3) 

where we are using Skyrme units in which lengths are measured in units of 2/{ef-^') and 
energy is measured in units of / 7 r/( 2 e). £4 is the Skyrme term and £0 is the baryon current 
density squared, which is inspired by the BPS Skyrme model [16]. The symmetry of the 
Lagrangian for 1/ = 0 is (5 = SU(2)lxSU(2)r acting on U ds U ^ U' = In the 

vacuum, this symmetry is spontaneously broken down to H SU(2)l+r, which in turn acts 
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on U as U ^ U' = gUg\ giving rise to the target space G/H ~ SU(2)l-r. The conventional 
potential term in the Skyrme model, viz. the pion mass term, is 1/ = m^tr ( 2 I 2 — U — W), 
and it explicitly breaks the symmetry G to SU(2 )l+r. 

In this paper, it will prove convenient to use the following notation of an 0(4) nonlinear 
sigma model where we express the field U in terms of four real scalar fields nA{x) {A = 
1, 2,3,4) with the constraint n\ = 1: 


U = inaCF°‘ + n 4 l 2 = n • t, 


( 4 ) 


where a = 1,2,3 is summed over, cr“ are the Pauli matrices and U‘^U = I 2 is equivalent to 
n • n = 1. We thus obtain the 0(4) sigma model with the Skyrme and six-derivative terms 


C = 
Ca = 


+ c^Ca -h cqCq - P(n), 

[(9^11 • d'^nf - (5^n • d^nf] , 

^ {e^^^^e^''P''nAdunBdpncd„nDf 

— ^ [(^/^n • cl^n)^ — 3 {dpO. • d^n) {dvii • dpuf -I- 2 {dpO. • d’'n) (d^n ■ 5^n) {dpU ■ d^n 


( 5 ) 

( 6 ) 
( 7 ) 
)]• 


The symmetry SO(4) ~ SU(2) x SU(2) for 1/ = 0 is thus manifest. 

The target space (the vacuum manifold with V — A) Ai — SU(2) ~ has a nontrivial 
homotopy group 


vr3(M) = Z, (8) 

which admits Skyrmions. The baryon number (the Skyrme charge) of S G 7r3(S'^) is defined 
as 

B = j {U^diUU^djUU^dkU) 

= ^ y {U^diUdjU^dkU) 

= j c^^^^e'-^'^dinAdjnBdkncnD 

= J d^x dinAd2nBd-incnD. (9) 

Instead of the conventional potential term, we consider here the following potential term 

V = rin?n\. (10) 
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We call this potential the Heisenberg type. The vacua of this potential is determined by 
724 = 0 , and thus the vacuum manifold is a sphere 


M-S^, ( 11 ) 

parametrized by Ui, n 2 , 72,3 with a constraint X]a=i 23 '^a ~ From the homotopy group 

7r2(M) ~ Z, (12) 

this model admits a monopole, viz. a global monopole. 

In our previous papers dH uni EDM, we considered instead the potential term V = 
m?(l — nl), i.e. a modified mass term. This admits two discrete vacua and a domain wall 
interpolating between them. 


In this paper, the potential (10) differs seemingly only by the overall sign. The physics 


described by this potential is, however, vastly different. The boundary condition at spatial 
infinity is 


77.4 = 0 , for r ^ 00 , 


( 13 ) 


which has the vacuum manifold E3 The vacuum state is thus any point on the two-sphere, 
given by 72^-1-722+? 2 | = 1. The potential breaks the symmetry G = SU( 2 )lxSU( 2 )r down 
to SU( 2 )L+i^, explicitly. This SU( 2 ) or 0(3) symmetry is, however, spontaneously broken 
down to 0(2). The spontaneous breaking 0(3)/0(2) gives rise to 2 Nambu-Goldstone bosons 
which remain massless as well as 1 massive pion. Abusing notation a bit, we can say that 
there are 2 massless pions and 1 massive pion. The phase of QCD we are trying to mimic, is 
high-density QCD at sufficiently high density; above the critical density for the formation of 


half-Skyrmions [H] and below the chiral restoration density. The potential (10) is an effective 
potential that can model this phase in the sense that it allows for half Skyrmions. It may 
or may not be the complete potential describing real QCD at mentioned densities, but in 
this paper we will study the above presented Skyrme-like models in the presence of said 
potential. 


III. FRACTIONAL SKYRMION MOLECULES 

In this section we will calculate numerical solutions of fractional molecules in the B = 1 
sector, which are bound states of two half baryons. The only free parameter is the mass ?72, 
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i.e. the coefficient in front of the Heisenberg-type potential. We will carry out the study in 
both the 2+4 and the 2-1-6 models. Instead of varying the mass, we will keep the mass fixed 
and vary the parameters (c 2 ,C 4 ) and (c 2 ,C 6 ) for the 2-|-4 and the 2-1-6 models, respectively. 
This includes the possibility of having C 2 = 0 which is not possible in the rescaled system 
where only the mass is varied. The C 2 = 0 region of parameter space is especially interesting 
due to the BPS properties that are present in that limit, see [16] for the 2-1-6 model and |33| 
for a possibility in the 2-|-4 model (although this requires a particular potential). 

As an Ansatz for the initial guess fed to the relaxation procedure, we will simply use the 
hedgehog Ansatz, suitable for the B = 1 sector 

n = (—cos/(r),xsin/(r)), (14) 

where x is the 3-dimensional spatial unit vector and = x^+y^ + z^ is the radial coordinate. 
The relaxation method will then deform the initial guess away from the spherical initial 
guess to the correct minimal-energy state of molecular shape (due to the Heisenberg-type 
potential). Note that we have chosen a particular value on the vacuum manifold, namely 
Hi = —1, which of course by 0(3) symmetry is equivalent to any other choice. 

We are now ready to perform the numerical calculation and we use the finite difference 
method to discretize the fields n, feed the initial guess to the relaxation algorithm and 
simply evolve a linear time operator until the numerical precision of the solution satisfies 
our criteria. 

We consider first the 2-1-4 model, i.e. C 2 > 0, C 4 > 0 and cq = 0. In Fig. [^is shown an array 
of molecule baryon charge density isosurfaces for various choices of the coefficients (02,04) 
for fixed mass m = 4. The coloring scheme used is based on the hue-saturation-lightness 
(HSL) parameters and the hue is given by the phase arg( 7 r 2 -|- ivrs), the lightness is given 
by IttiI and 7 ? = n/\n\ is a normalized 3-vector. All the numerical calculations throughout 
the paper are carried out on an 81^ cubic lattice using the relaxation method. We observe 
that the molecular shape is most pronounced when the coefficients are small (which is clear, 
because that corresponds to a large mass) and even more so when 02 04. 

In Figs. 1^ and are shown cross sections at .2 = 0 of the baryon charge density and 
energy density, respectively. The numerically integrated baryon charge density, denoted by 
^numerical g^ygg ^ handle on the precision of the numerical solution, see Tab. |l| Note that 
the molecular shape is slightly more pronounced in the energy density than in the baryon 


7 




C2 = 0, C4 = I 


C2 = 0, C4 = 1 


C 2 = 0, C 4 = 4 








C 2 = 1, C4 = I C 2 = 1, C4 = 1 C 2 = 1, C 4 = 4 

FIG. 1: Isosurfaces showing the half-maximum of the baryon charge density in the 2-h4 model 
for various choices of ( 02 , 04 ) for fixed mass m — A. The color scheme represents the normalized 
3-vector tt = where arg( 7 r 2 + ins) is the hue and the lightness is given by | 7 ri|. 



charge density, viz. the depth of the valley between the two peaks is deeper. In the cases 
02 = 1 , 04 = I and 02 = 04 = 1 , there is no valley between the two peaks; there is however 
still some amount of “polarization;” or better, there is a dipole moment. In order to quantify 










the amount to which the Skyrmion is moleculized, we define the following quantity 


= 2 d^x (sign(n4)xi5 — sign{n3)yB) , 


( 15 ) 


which we call the baryonic dipole moment and 

1 


B = 


^ABCD 


2'k‘^ 


dinAd2nBdincnD, 


(16) 


is the baryon charge density. The physical meaning of which has units of length, is 
to which degree the baryon charge corresponding to the northern and southern hemisphere 
(distinguished by the sign of n^) is separated compared to how the charge is distributed in 
the transverse direction. This quantity is absolute and should be compared to the size of the 


Skyrmion (see the figures). For a Skyrmion without the potential (10), the baryonic dipole 
moment vanishes: p^ = 0. For convenience we define the size of the Skyrmion as 
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B _ 


(Px 


(17) 


We mentioned already that the two coefficients 02,04 are not independent of each other 
and they can be scaled to unity giving a mass which we will denote TT^canomcai can be 

calculated as 


^canonical 



C2 


(18) 


where m is the mass in the noncanonically normalized Lagrangian, i.e. with C 2 7 ^ 1 and 
C 4 7 ^ 1. The values of for the various configurations shown in Figs. 01 and 1 ^ are 

given in Tab. |Ij 

We will now consider the 2+6 model, i.e. 02 > 0, Oe > 0 and 04 = 0. In Fig. is shown 
an array of molecule baryon charge density isosurfaces for similar values of the coefficients 
( 02 , 00 ) as in the previous case and for fixed mass m = 4. The coloring scheme is the same 
as used above. Again the molecular shape is mostly pronounced when the coefficients are 
small and 02 -C og. We also calculate the canonical mass, 


^canonical _ 


i.e. the mass corresponding to a rescaled equation of motion with 02 
those corresponding values in Tab. |TTj 


= Cg = 


(19) 

1 and show 
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C 2 = 1, C 4 = 1 C 2 = 1, C 4 = 4 

FIG. 2 : Baryon charge density at a spatial slice through the molecule at 2 ; = 0 in the 2+4 model 
for various choices of ( 02 , 04 ) and for fixed mass m = 4. 



In Figs. and 1^ are shown cross sections at 2 ; = 0 of the baryon charge density and energy 
density, respectively. The nnmerically integrated baryon charge density, denoted ^“umericai 
gives again a handle on the precision of the numerical solution, see Tab. [TT| We observe that 
the energy density is more spiky than the baryon charge density, hence more of a molecular 
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25 


25 



C2 = 1,C4 = ^ C 2 = 1, C 4 = 1 

FIG. 3: Energy density at a spatial slice through the molecule at z 
various choices of ( 02 , 04 ) and for fixed mass m = 4. 


C 2 = 1,C4 = 4 
0 in the 2+4 model for 


shape. 

The case of C 2 = 0 and Ce > 0 is particularly interesting because of its BPS and integrable 
properties. It is by now called the BPS Skyrme model [16]. The integrable property is very 
appealing since analytic solutions can readily be calculated for a large class of potentials. 
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TABLE I: The numerically integrated baryon charge, the numerically integrated energy and the 
numerically integrated baryonic dipole moment for the various configurations in the 2+4 model, 
shown in Figs. [Il[2]and|3| 


C2 

C4 

^canonical 

^numerical 

^numerical 


0 

1 

4 

OC 

0.99969 

18.930 

0.575 0.772 

0 

1 

oc 

0.99986 

52.720 

0.970 1.151 

0 

4 

OC 

0.99989 

154.74 

0.998 1.472 

1 

4 

1 

4 

23 

0.99989 

25.832 

0.288 0.620 

1 

4 

1 

2^ 

0.99967 

65.527 

0.478 0.901 

1 

4 

4 

23 

0.99989 

167.24 

0.920 1.437 

1 

1 

4 

2^ 

0.99962 

40.078 

0.087 0.426 

1 

1 

22 

0.99958 

89.365 

0.185 0.713 

1 

4 

23 

0.99972 

203.93 

0.635 1.282 


TABLE II: The numerically integrated baryon charge, the numerically integrated energy and the 
numerically integrated baryonic dipole moment for the various configurations in the 2+6 model, 
shown in Figs. m [5] and m 


C2 

C6 

^canonical 

^numerical 

^numerical 


1 

4 

1 

4 

23 

0.99989 

34.440 

0.595 0.978 

1 

4 

1 

2i 

0.99986 

61.728 

0.863 1.278 

1 

4 

4 

2^ 

0.99986 

113.89 

1.353 1.707 

1 

1 

4 

2I 

0.99986 

60.705 

0.281 0.822 

1 

1 

22 

0.99984 

97.104 

0.440 1.091 

1 

4 

2I 

0.99985 

161.15 

0.708 1.442 


Numerically, however, it is a rather difficult problem, because the Skyrmion in the BPS limit 
turns into a compacton cni, giving the soliton a finite size (hence no exponentially damped 
tail) and so a cusp in the fields at a finite distance. There is no cusp in the energy density, 
which is smooth, but the cusp in the fields requires a special technique in order to be studied 
numerically. We will not pursue this problem further in this paper. 

The half-Skyrmion molecule has two components which we can interpret as global 
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C2 — 1, C6 — I C2 — 1, C6 — 1 C2 — 1, C6 — 4 


FIG. 4: Isosurfaces showing the half-maximum of the baryon charge density in the 2+6 model for 
various choices of (c 2 , cq) for fixed mass m = 4. The color scheme is the same as that in Fig. 

monopoles, both of baryon charge 1/2. In Fig. is shown a cross section at z = 0 of 
the 714 component of the molecule in the 2+4 model with C 2 = ^ and C 4 = 1, as an exam¬ 
ple. The same characteristic holds for all the molecules that we found in this section. The 
figure nicely shows that 77.4 = 0 , not only at infinity, but also on a plane separating the two 
constituents of the molecule (which is seen as a line in the cross section shown in Fig. [^. 
In the next section we will construct only one of the two constituents isolated, however, at 
the cost of having an infinite total energy of the configuration. 

IV. FRACTIONAL SKYRMIONS AS GLOBAL MONOPOLES 

In this section we will consider a limit in which one constituent of the molecule has been 
drawn away to infinity, thus leaving the molecule with only a half baryon number. The result 
is a half Skyrmion. This is usually not possible. The loop hole is that the configuration has 
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C2 — 1, C6 — I C2 — 1, C6 — 1 C2 — 1, C6 — 4 

FIG. 5: Baryon charge density at a spatial slice through the molecule at z = 0 in the 2+6 model 
for various choices of (c 2 , cq) and for fixed mass m = 4. 


a divergent total energy. By inserting the Ansatz 


n = (xsin/(r),cos/(r)), 


( 20 ) 


into the Lagrangian ([^, we get 


C2 


C2 


C4 


C4 


Cq 


m 


^ sin / + - sin (/)/^ + ^ (21) 


which has the corresponding equation of motion 

C2 f /rr + -/r - 4 siu 2/^ + ^4 f 4 Sm^(/)/rr + 4 ^M‘^f)fr “ 4 f 2/ 


rp p2 ^ I \ ' p2 \ j JJr ^4 

+•^6 sin^{f)frr - ^ sin^{f)fr + ^ sin^ / sm{2f)ff^ + sin 2/ = 0, (22) 

where = drf. The difference between this Lagrangian with corresponding equation of 
motion and the normal case with a mass term is that the potential sin^ / is replaced with 
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FIG. 6 : Energy density at a spatial slice through the molecule at z = 0 in the 2+6 model for 
various choices of (c 2 , ce) for fixed mass m = 4. 



FIG. 7: The 77.4 component at a spatial slice through a molecule at 2 ; = 0 in the 2+4 model with 
C 2 = I and C 4 = 1 . 
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cos^ / and in turn — sin 2/ in the equation of motion becomes +sin2/. We therefore need 
to consider the boundary conditions for this system in order to construct a half Skyrmion. 
One possibility is to choose 

TT 

southern molecule constituent: /(O) = tt, f{oo) = —, (23) 

giving a half Skyrmion winding on only the southern hemisphere of the target space. Alter¬ 
natively, we could have chosen 


northern molecule constituent: /(O) = 0, 


/(oo) = 


TT 

2 ’ 


(24) 


giving instead a half anti-Skyrmion which winds only on the northern hemisphere of the 
target space. 

In Figs. 1^ and are shown numerical solutions for the half Skyrmion winding on only 
the southern hemisphere of the target space (i.e. n 4 G [—1,0]) for the 2-|-4 and 2-1-6 model, 
respectively. The half Skyrmion can live in isolation only at the cost of an inhnite total 
energy. Figs. and show that the energy density multiplied by r^, i.e. goes to 

a constant and hence the total energy picks up a linear divergence in the radial integral 
E oc where i?max is the radial cut off of the integral. All the baryon charge densities 

integrate numerically to one half. 


V. HIGHER BARYON NUMBERS 


In this section we make a first attempt to make solutions with higher baryon numbers, 
i.e. B > 1. As we ab initio do not know the shape of the global minimizers, we use the 
relaxation method with different initial guesses. The two guesses we choose here are the 
axially symmetric Ansatz 


n = (—cos/(r), sin6*sinS0sin/(r), sin6*cos S0sin/(r), cos6*sin/(r)), 


(25) 


and the following Ansatz using a rational map R: 

( f( \ ■ f! \ . f( \ '^ — RR . , 

n = I — cos / (r),-—= sm j (r),-—sm / (r), ^ ^ - sm j 


1 + RR 


1 + RR 


1 + RR 


(26) 


with the symmetries found to minimize the Skyrmions without the potential (10). For the 


appropriate rational maps, R, see Refs. 1201 ED- Note again that we have chosen a particular 
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FIG. 8: An isolated half Skyrmion (a global monopole) in the 2+4 model for varions masses m: 
(a) the radial field profile = cos/, (b) and (c) the energy density S and 47rr^£^, respectively. 
The latter shows the divergence in the total energy, (d) is the baryon charge density. The chosen 
parameters are C 2 = \, C 4 — 1 and m — 1,2,3,4. 


value on the vacuum manifold, i.e. ni = —1 which by 0(3) symmetry is equivalent to any 
other choice. 

We begin with the 2+4 model and calculate the numerical solutions for the first six 


baryon numbers, i.e. 5 = 2,3,4, 5 ,6 (5 = 1 was already made in Sec. III). For concreteness, 
we fix the parameters C 2 = |, C 4 = 1 and m = 4, corresponding to a canonical mass 


m 


canonical _ 


= 16. We find that the lowest-energy states take the shapes of rings (beads on 
rings) and show the numerical solutions as isosurfaces at the half-maximum baryon charge 


densities in Fig. 10. These solutions are made with the axially-symmetric initial guess (25). 
We again use the numerically integrated baryon charge density, ^numerical handle on the 
precision of the numerical solution, see Tab. ill In this table we also show the total energy 
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FIG. 9: An isolated half Skyrmion (a global monopole) in the 2+6 model for varions masses m: 
(a) the radial field profile = cos/, (b) and (c) the energy density S and 47rr^£^, respectively. 
The latter shows the divergence in the total energy, (d) is the baryon charge density. The chosen 
parameters are C 2 = \, ce — I and m — 1,2,3,4. 


per unit baryon charge, ^’numencaiy'^^ which tells us whether the higher-charged solution is 
stable or only metastable. We find that all the higher-charged solutions we have calculated, 
namely B = 2,3, 4, 5,6 are in fact stable (at least among this type of configurations). 

The isosurface shows the three-dimensional structure of the solution, but not the profile 
shape of the baryon-charge density or energy density. Therefore, as before, we show cross 
sections at z = 0 of the baryon charge density and energy density in Fig. [T^ and 
respectively. Note that again the molecular shape is slightly more pronounced in the energy 
density than in the baryon charge density, viz. the depth of the valleys between the peaks 
are deeper. Finally, we notice that the molecular shape (by which we mean that half a unit 
of baryon charge is spatially localized) is far more pronounced for higher baryon charges 
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FIG. 10: Isosurfaces showing the half-maximum of the baryon charge density in the 2+4 model for 
baryon numbers B = 2,3,4, 5,6 with C 2 = |, C 4 = 1 and fixed mass m = 4. The color scheme is 
the same as that in Fig. [TJ 

B > 1 than for S = 1, for the same coefficients, i.e. C 2 = C 4 = 1 and m = 4. 
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TABLE III: The numerically integrated baryon charge and energy per unit baryon charge for higher 
baryon numbers in the 2+4 model. 


^ ^numerical ^numerical j^ 


0.99967 

65.52(7) 

1.9994 

61.5(3) 

2.9979 

61.0(6) 

3.9971 

60.8(9) 

4.9964 

60.7(9) 

5.9957 

60.5(5) 


As we mentioned already, we do not ab initio know the symmetries or spatial structure of 
the global energy-minimizing solutions. Therefore we have tried also different initial condi¬ 
tions to check whether we can obtain lower-energy solutions with the same baryon numbers 


compared to those obtained with the axially-symmetric initial guess (25). Explicitly, we 
have tried the rational-map Ansatz with the rational maps that minimize the energy with¬ 
out the potential iH). We found these solutions to have higher energies than those shown 
here (or being spatially disconnected) and more details are given in Appendix |Aj 

We will now consider the 2-1-6 model and calculate numerical solutions for the hrst five 


baryon numbers, i.e. S = 2,3,4, 5 (5 = 1 was made in Sec. III). For concreteness, we fix here 
the parameters C 2 = |, Ce = 1 and m = 4, corresponding to a canonical mass jTicanonicai _ 
8-^/2. As in the case of the 2-1-4 model, we find again that the lowest-energy states take 
the shapes of rings (beads on rings). These solutions are made with the axially-symmetric 


initial guess (25). In Fig. 13 are shown the isosurfaces of the baryon charge densities at their 
respective half-maximum values of the numerical solutions. The numerically integrated 
baryon charge and total energy per unit baryon charge are displayed in Tab. [IVl We find as 
in the case of the 2-1-4 model, that all the calculated solutions with higher-baryon numbers, 
B > 1, are stable (as opposed to metastable) among this type of configurations. 

The cross sections at = 0 of the baryon charge densities and energy densities are 


shown in Fig. p^andfl^ respectively. Note that again the molecular shape is slightly more 
pronounced in the energy density than in the baryon charge density, viz. the depth of the 
valleys between the peaks are deeper. In fact, the slices of the energy densities suggest that 
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FIG. 11: Baryon charge density at a spatial slice through the B — 2,3,4, 5, 6 molecules at z = 0 
in the 2+4 model for various choices of (c 2 , C 4 ) and for fixed mass m — 4. 


the beads are not spatially connected. This is, however, not true. From the slices of baryon 


charge densities, we see that the beads are connected and from Tab. [TV] we can see that 
there is in fact a small binding energy. Asymptotically, however, for i? —>■ 00 the binding 
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FIG. 12: Energy density at a spatial slice throngh the B = 2,3,4, 5, 6 molecules at z = 0 in the 
2+4 model for various choices of ( 02 , 04 ) and for fixed mass m — 4. 


energy may go to zero (in the 2+6 model for the given parameters). 
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FIG. 13: Isosurfaces showing the half-maximum of the baryon charge density in the 2+6 model for 
baryon numbers B = 2,3,4,5 with C 2 = \, cq = 1 and fixed mass m = 4. The color scheme is the 
same as that in Fig. [T] 


VI. MOLECULES WITH UNEQUAL FRACTIONS 


In this section, we show that by a modification of the potential (10), we can create a 
molecule, still with two components, but with unevenly distributed baryon charge. The 
modified potential is 




(27) 


giving rise to a molecule with two components with localized baryon charge (1 + c)/2 and 
(1 —c)/2 , respectively. This potential thus requires different asymptotic boundary conditions 
and hence different initial guesses. Since we consider only the B = 1 sector here, a modified 
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TABLE IV: The numerically integrated baryon charge and energy per unit baryon charge for higher 
baryon numbers in the 2+6 model. 


^ ^numerical ^numerical j^ 


1 0.99986 61.72(8) 

2 1.9998 59.2(7) 

3 2.9990 59.1(6) 

4 3.9985 59.0(7) 

5 4.9981 58.9(1) 






FIG. 14: Baryon charge density at a spatial slice through the B = 2,3,4,5 molecules at z = 0 in 
the 2+6 model for various choices of (c 2 , ce) and for fixed mass m — A. 
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FIG. 15: Energy density at a spatial slice through the B = 2,3,4, 5 molecules at z = 0 in the 2+6 
model for various choices of (c 2 , ce) and for fixed mass m = 4. 


hedgehog Ansatz is sufficient 

n = sin /(r) — \/l — cos f {r),y sin f{r), zsin/(r), ccos /(r) + \/l — c^x sin /(r)^ , 

(28) 

which satisfies the vacuum equation at r oo provided f{oo) = 0 . 

In the following we will choose a concrete example, setting c = |, which is interesting, 
because the two components of the baryon should contain | and | of the unit baryon charge, 
respectively For concreteness, we will keep the parameters that we have used in Sec. |V| 
namely: C 2 = |, C 4 = 1 (ce = 1) and m = 4 in the case of the 2+4 (2+6) model. The 


numerical solutions are shown in Figs. 16 and 17 
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FIG. 16: Unequal fractional molecule with c = ^ giving a molecule in the 2+4 model with two 
components of charge | and respectively, (a) shows the isosurface of the baryon charge density 
at half-maximum value and there is an added shadow which is an isosurface at one quarter of the 
maximum value, (b) and (c) show x^-slices (at z = 0) of the baryon charge density and the energy 
density, respectively. The parameters are C2 = C4 = 1 and m = 4. The color scheme is the same 
as that in Fig. The numerically integrated baryon charge is ^numerical _ 0.99974. 



(a) isosurface (b) baryon charge density (c) energy density 

FIG. 17: Unequal fractional molecule with c = ^ giving a molecule in the 2+6 model with two 
components of charge | and respectively, (a) shows the isosurface of the baryon charge density 
at half-maximum value and there is an added shadow which is an isosurface at one quarter of the 
maximum value, (b) and (c) show xy-slices (at 2; = 0) of the baryon charge density and the energy 
density, respectively. The parameters are C2 = |, ce = 1 and m = 4. The color scheme is the same 
as that in Fig. The numerically integrated baryon charge is ^numerical _ 0.99995. 
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VII. SUMMARY AND DISCUSSION 


We have constructed fractional Skyrmions and their molecules in the Skyrme model 
with the potential term, V = m?n^. As for higher-derivative terms, we have considered the 
conventional fourth-order derivative term, viz. the Skyrme term or the sixth-order derivative 
term being the baryon number current squared. One molecule consists of a pair of a global 
monopole and anti-monopole, each with half a baryon number. Since an isolated global 
monopole has divergent energy in an infinite space, the constituents are confined. We have 
also constructed an isolated fractional Skyrmion as a global monopole (which thus has a 
divergent total energy). We have then constructed Skyrmion solutions with higher baryon 
numbers up to 5 = 6, and have found that configurations in the form of beads on rings 
are energetically stable. We have also found other metastable configurations for S = 3 and 
B = 5, but exhausting all possible metastable configurations remains as a future problem. 
Finally by considering the potential term V ~ m?{n^ ~ we have found that fractional 
Skyrmions have baryon numbers that are not equal to one half. As an example, we have 
constructed a molecule with fractional Skyrmions with the baryon numbers 1/3 -|- 2/3. 


We have shown that our choice of the potential (10) as an effective low-energy potential is 
able to describe the half-Skyrmion phase which should be present in QCD at high density [6]. 
In Ref. [0], which studies a holographic model, quite different from our low-energy effective 
field theory model, the u) and p mesons were needed to consistently describe the high density 
phase where the half Skyrmions exist. Whether our simple model with our choice of effective 
potential is able to capture the relevant phenomenological features of QCD at high density 
is a very interesting and important question that is left as a future work. 

As a lower dimensional analog, there exists a fractional baby-Skyrmion molecule consist¬ 
ing of a pair of a global vortex and anti-vortex with half 1^2 charges in the 0(3) model with 
the XY (or easy-plane) potential term V = in d = 2 -|- 1 dimensions. In this model 

too, baby Skyrmions as beads on a ring were found as (meta)stable configurations [TUI ITT] . 
However, the lowest-energy configurations of higher topological numbers are of the form of 
square lattices of fractional lumps [9]. We expect a cubic lattice of fractional Skyrmions for 
higher baryon numbers in our 3-1-1 dimensional case. 

If we add a potential V 2 = nri^ni with m 2 < m in the 0(3) model with the XY-potential 
term V = in d = 2 -|- 1 dimensions, vortices are connected by a sine-Gordon soliton. 
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Therefore, fractional lumps are linearly confined by the sine-Gordon soliton. This happens 
in fact in two-gap (or two-component) superconductors described by a Landau-Ginzburg 
Lagrangian in the form of a U(l) gauge theory with two charged scalar fields 01 and 02- 
The two fields are coupled through the Josephson term Lj = 70102 in the case of two-gap 
superconductors, but not for two-component superconductors. In either case, it admits a 
molecule of half-quantized vortices [3l]. In the presence of the Josephson term, fractional 
vortices are connected |3S] by a sine-Gordon kink [SH]- In the strong gauge-coupling limit, 
the model reduces to an 0(3) sigma model with the potential term V = complemented 
by 14 = yni induced by the Josephson term (33 Ei, and the molecule reduces to a fractional 
baby-Skyrmion molecule mentioned above [39] . In the same way, if we add a potential that 
breaks the S0(3) symmetry, possessed by the vacuum, such as 14 = m^ni (m 2 < m), 
fractional Skyrmions constituting a molecule will be connected by a baby-Skyrmion string, 
realizing linear confinement of fractional Skyrmions. 

The Bogomol’nyi-Prasad-Sommerfield (BPS) Skyrme model, proposed recently (Hj , con¬ 
sists of only the sixth-order derivative term as well as appropriate potentials. This model 
admits exact solutions with compact support. By choosing the potential in this paper, we 
may be able to construct exact solutions of a fractional Skyrmion molecule. 

Fractional Skyrmions in the 0(4) model (or the Skyrme model) on x were discussed 
in Ref. (30] , in which our potential term is related to the boundary condition where the field 
n 4 changes the sign along S^. A rather different origin of fractional topological charge was 
also found for vortices and lumps [H]. A unified understanding of fractional topological 
charges will be an important future work. 

Acknowledgments 

The work of M. N. is supported in part by Grant-in-Aid for Scientific Research 
No. 25400268 and by the “Topological Quantum Phenomena” Grant-in-Aid for Scientific 
Research on Innovative Areas (No. 25103720) from the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT) of Japan. S. B. G. thanks Keio University for 
hospitality during which this project took shape. S. B. G. thanks the Recruitment Program 
of High-end Eoreign Experts for support. 


28 


Appendix A: Rational map Ansatz in the 2+4 model 


In this section we investigate numerical solutions with different initial guesses to be used 


by the relaxation method, in particular, we consider the rational-map Ansatz (26) with an 


appropriate rational map R. For B = 1 and B = 2, the axially symmetric Ansatz (25) is a 
very educated guess due to the high level of symmetry. For B = 3 and higher B, without 


the potential ( 10 ), the lowest-energy state turns out to have a discrete symmetry instead of 
axial symmetry 1201121]. As already mentioned, we do not a priori know what symmetry 
or shapes the minimizer of the energy may have, and therefore we try the minimizers which 


are found to be the lowest-energy states without the potential ( 10 ) as initial guess also here. 


i.e. with the potential (10). For concreteness, we use the 2+4 model with C 2 = ^, C 4 = 1 
and m = 4 and calculate numerical solutions for = 3 ,4, 5. The isosurfaces of the baryon 
charge density at half-maximum values are shown in Fig. 

The first case, B = 3, found using the rational map with a tetrahedral symmetry turns 
out to give a numerical solution which is only metastable, see Tab. |V| for a comparison of 
the energies of the numerical solutions. In order to demonstrate that the two clusters of 
three half units of Skyrme charge are really connected, we display a cross section of the 
baryon charge density and energy density at a: = 0 in Fig. Of course, no half unit of 
Skyrme charge can be spatially localized (without a tail) and thus neither can three half 
units. Both the baryon charge density and the energy density on the yz slices in Fig. [T^ 
have maximum values around one quarter of their respective global maximum values (i.e. in 
the three-dimensional space). 

The case, 5 = 4, found using the rational map with a cubic symmetry gave rise to 
a solution that split up into two 5 = 2 axially symmetric solutions; i.e. they are both 
chargewise and energetically separated. Due to the 5 = 2 solution with axial symmetry 
having a larger energy per unit baryon charge than the 5 = 4 solution, see Tab. in two 
5 = 2 solutions thus have a higher energy than the 5 = 4 solution being a ring (with 
beads). 

The last case, i.e. 5 = 5 is like the 5 = 3 case spatially connected, but is made with a 
rational map having an octahedral symmetry. It is however energetically only metastable, 
viz. it has a higher energy than the solution made using the axially symmetric Ansatz. 

The bottom line is that the lowest-energy solutions for the molecules are all found using 
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FIG. 18: Isosurfaces showing the half-maximum of the baryon charge density in the 2+4 model for 
baryon numbers B = 3,4, 5 with C 2 = ^, 04 = 1 and fixed mass m = 4. These numerical solutions 


have been obtained using the rational-map Ansatz (26) as initial guesses. The color scheme is the 
same as that in Fig. 


the axially symmetric Ansatz (25). Using the rational-map Ansatz (26), we have found 
metastable solutions for the B = 3,5 cases, which however are energetically prone to decay. 


Appendix B: Low-mass limit of the 2+4 model 

For completeness, we provide a series of numerical solutions for the 2+4 model, for the, in 
this paper, most studied case; i.e. C 2 = |, C 4 = 1 with various masses m = 0, 1, |, 2, |, 3, | 

whereas the case of m = 4 is used throughout the paper. The isosurfaces of the baryon 
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TABLE V: A comparison of the numerically integrated baryon charge and energy per unit baryon 
charge for higher baryon numbers in the 2+4 model for different initial guesses. The denotes a 

solution that is spatially disconnected. 

B (initial guess) ^numerical £;numericay^ 


3 (axially symmetric) 2.9979 

61.0(6) 

3 (tetrahedral) 

2.9986 

61.1(7) 

4 (axially symmetric) 3.9971 

60.8(9) 

4^ (cubic) 

3.9971 

61.5(2) 

5 (axial) 

4.9964 

60.7(9) 

5 (octahedral) 

4.9964 

61.5(4) 




baryon charge density energy density 

FIG. 19: Cross sections at x = 0 showing the baryon charge density and energy density in the 
2+4 model for baryon numbers B = 3 with C 2 = |, C 4 = 1 and fixed mass m — A. This numerical 


solution has been obtained using the rational-map Ansatz (26) as an initial guess. 


charge densities of the solutions are shown in Fig. The baryon charges, masses, dipole 


moments and sizes are given in Tab. VI 
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m = 2 m = ^ m = 3 m = ^ 

FIG. 20: Isosurfaces showing the half-maximum of the baryon charge density in the 2-h4 model for 
a single baryon with masses m = 0, 1, 2, 3, ^ and the usual fixed coefficients C 2 = |, C 4 = 1. 

The color scheme is the same as that in Fig. [T] 


TABLE VI: The numerically integrated baryon charge, the numerically integrated energy and the 


numerically integrated baryonic dipole moment for the various masses shown in Fig. 20 


m 

^numerical ^numerical 



0 

0.9963 

35.731 

0 

1.876 

1 

2 

0.9895 

38.044 

0.334 

1.763 

1 

0.9997 

43.930 

0.364 

1.441 

3 

2 

0.9999 

48.457 

0.398 

1.271 

2 

0.9999 

52.446 

0.419 

1.157 

5 

2 

0.9998 

56.069 

0.432 

1.072 

3 

0.9998 

59.421 

0.439 

1.005 

7 

2 

0.9997 

62.556 

0.443 

0.950 

4 

0.9997 

65.527 

0.478 

0.901 
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